Markovian evolution of Gaussian states in the semiclassical limit 
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We derive an approximate Gaussian solution of the Lindblad equation in the semiclassical limit, 
given a general Hamiltonian and linear coupling with the environment. The theory is carried out 
in the chord representation and describes the evolved quantum characteristic function, which gives 
direct access to the Wigner function and the position representation of the density operator by 
Fourier transforms. The propagation is based on a system of non-linear equations taking place in 
Q\ , a double phase space, which coincides with Heller's theory of unitary evolution of Gaussian wave 

packets when the Lindbladian part is zero. 

■ 

PsJ ' PACS numbers: 03.65.Sq,03.65.Yz 

; I. INTRODUCTION 

: 

Open quantum systems deal with physical systems which interact with their environment. From this interaction 
the system may loose or gain energy; this is the dissipation phenomenon. Coupling of a system with its environment 
entails loss of information in the sense that, if one starts with a pure state p(0), such that its purity, Trp(0) 2 = 1, 
then it will undergo a non-unitary evolution which will not preserve its purity, that is, Trp(t) 2 < 1. One important 
aspect of this loss is decoherence, that is the vanishing of the off-diagonal terms of the density operator. 

The Markovian process offers a powerfull mathematical tool in the study of this class of systems, ft relies on 
the assumption that the environmental degrees of freedom are very fast as compared to the proper dynamics of the 
system, so that the future evolution of the density operator is a function of its value in the present only, regardless of 

1 ^ i ' 1 its past. It was shown by a series of works, concluded by Lindblad in that the corresponding evolution equation 

can always be written in the form 

m ■ ^ 

in ! k 

The Hamiltonian H describes the unitary evolution of the system without environment, and Lindblad's theory proves 
qq , the a priori existence of the operators Lk, which are commonly dubbed Lindblad operator, and which modelize the 
effect of the environment. The master equation of quantum optics is a particular, well understood case, where the 
Lindblad operators are L\ — a, i.e. the annihilation operator, describing the emission process, and L2 = cr the 
absorption. One can see that these Lindblad operators are not Hermitian, which can be generalized to every case 
where the coupling to the environment is dissipative. 

In this paper, we derive the consistent dynamics of a Gaussian solution of the Lindblad equation. As compared to 
former paper where we generalized the analysis of the evolution of extended states, that goes back to Van Vleck \m 
(see also [7] 8} 9\), this work is rather a generalization of the evolution of "wave packets", developed by Heller [l(| 
and Littlejohn [11|, among others. The approximation holds as long as the size of the Gaussian is small enough to 
justify an identification of the Hamiltonian with its second order expansion, that is, its quadratic kernel. Since the 
effect of decoherence is to quickly decrease the extension of the solution, the regime of validity should be larger than 
in the unitary case. 

We chose to represent the Gaussian solution in the chord reresentation, that is, the Fourier transform of the Weyl- 
Wigner representation. After giving the form of the Lindblad equation in the chord frame, we assume a Gaussian 
expression of the solution, with time dependent parameters, and its insertion in the Lindblad equation leads to a 
system of equations of motion for these parameters. Comparison with Heller's theory is then derived in a slightly 
simplified case. 

This paper presents a simpler, albeit less exact theory than [ijj], which generalizes complex WKB approximations 
for the solutions of the Lindblad equation. In that treatment it is necessary to further complexify the doubled phase 
space. Thus, it is likely that the simpler approximations developed here may be of more practical use. 

The theory can be useful for any initial state, which can be conceived as a finite decomposition into Gaussian ones. 
All the following formulae are written for a system with a single degree of freedom, in order to clarify the notation. 
Nonetheless, it is quite simple to generalize them for a finite number of degrees of freedom. 
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II. EQUATION IN THE WEYL FORMALISM 



The Weyl representation maps every quantum operator onto a phase space function, that is, a function of the vector 
x = (P: l) [IB] [13] • For an operator A the Weyl symbol A is defined as 



Mx)=2 j exp(-^Q) (q+Q\A\q- 



The Weyl symbol of the state operator p is the Wigner function 



-} dQ. 



W(x)=N / exp(--pQ) (« + ^|p|«--^) dQ, 



with A/" = l/(27r7i), and its Fourier transform, the chord function x(£)i a ^ so called characteristic function, is 

X(0 = M I exp (--£ A cc) W(a?) ete, 



(2) 



(3) 



(4) 



where the wedge product of two vectors x = (j>, q) and x' — (p 1 , q') is defined by x A x' = pq' — p'q = Jx ■ x', which 

1 One can have the chord symbol directly from the quantum operator 



also defines the skew matrix JJ = 
through the formula 
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\(£>=A'/<'M>(™&9) (<? + fl/%-f> «/./• 



(5) 



We call the space of all x the centre space, and the space of all £ the chord space. 

In the chord space, by using product rules for the product of operators [17|, the Lindblad equation is represented 
by a partial differential equation. This equation is actually simpler than in the Weyl (centre) representation, and this 
justifies our choice. In the case where the Lindblad operators are linear functions of p and q, that is, L = I ■ x + il" ■ x 
with I' and I real vectors, this equation can be written 



= ~tf[[H& + ~e,t)-H&-~t,t)] exp A*') x (£'./),/£' ,/*' 



The dissipation coefficient, 



7 = Z"AZ', 



(6) 



(7) 



is null for a Hermitian Lindblad operator (Z = 0) and we then have a purely diffusive case. H is the Weyl represen- 
tation of the Hamiltonian of the isolated system and coincides with the corresponding classical Hamiltonian, up to 
corrections coming from non-commutativity of p and q. Its arguments in equation ^ are the pair of remarkable points 

x + = x+l and x =x- i which can be considered as both tips of a chord €; Although this chord € = ( tp , Q can 
be interpreted as an auxiliary conjugate variable of x, in the current approach it is actually more convenient to write 
the solution in terms of y = J£ = (— £ q ,£ p ). The direct sum of these conjugate spaces can indeed be interpreted as 
a double phase space, where x formally plays the role of the position g, and y the role of its Fourier conjugate, the 
momentum p. Then the above equation becomes 



^(y, t) = - % -N j [H(x' - ijy, t) - H(x' + ijy, t)] exp - y) ■ x' ) X ( y 1 . 1 ! dy da' 

-7 y ^(<W) 2 x(y,t). 



(8) 



The same name has been kept for the characteristic function %(y, i), though strictly this should be i) = Jl/j 
and we have set A = J(Z' + iZ"). The term t/- ^ can actually be included in the integral term, thanks to an integration 
by parts of the exponential, and one has finally 

^{y,t) = - l -N f H(x',y,t) expf-Cy'-y)-*') X t{v') dy 1 dx' i(A • yf X (v,t), (9) 
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with 



H{x', y, t) = H(x' - ijy, t) - H(x' + ijy, t) - 1 x' ■ y 
= H + (x',y,t)-H-(x',y,t)-f x' y. 



(10) 



This is exactly the double Hamiltonian that generates the classical motion underlying the semiclassical approximations 
in [IH . Obviously, the double hamiltonian will be time- independent if it is obtained from a time- independent single 
Hamiltonian. Furthermore, in the absence of dissipation, both Tt + (x, y) = H(x — \H>y) and Tt~(x, y) = H(x + ^Jy) 
will also be constants which generate independent motions for both chord tips. 
Notice that ([9]) can also be written as 



^(y,t) = -~N J H(x',y,t) exp(-~y ■ x') W t (x>) dx' - i(A- yf x (v,t), 



in terms of the evolving Wigner function, or, alternatively, as 



%t(v,t) = -^w(-?# (1) ,y (2) ,t) x(y,t) - ^(A-y) 2 X (y,t), 



dt 



i dy 



2h" 



(11) 



(12) 



where W and ^ mean that the derivatives are taken first and then the y terms are multiplied. 

The differential term in the RHS of (fT2|) (or the integral term in the RHS of (fTTj) ) represents the unitary part of 
the evolution. In other words, in the chord representation, the commutator is 



H,p 



= H{--^- ,yW,t) x (y,t). 

chord I Oy 



(13) 



III. GENERAL DYNAMICS 



We here assume that the chord representation x(y, t) of the localized state has the usual semiclassical form 

\{y,t) = exp-S(y,t), 



(14) 



where S(y,t) is a function with complex values of order 0(h ). It is important to notice that, if we find the time 
evolution of such a state determined by the Lindblad equation, then we can evolve any linear combination of such 
states, which can be coherent states for instance. This is a consequence of the linearity of the Lindblad equation. 

This semiclassical form naturally induces an h expansion of the unitary part of the equation, as it is shown in the 
appendix El 



H(--^- {1 \y^,t) X (y,t) 
i dy 



x(y,t). 



Hence, by expanding (|12[) at leading order in h, one obtains the Hamilton Jacobi equation 



OS, . 



-H 



^(y,t\y,?j+^(X-y) 2 + 0(h). 



(15) 



(16) 



This is a double phase space generalization of the complex WKB theory in [18[. In this paper, we will not develop 
a complex resolution of this equation, presented in [l9j]. We rather separate S(y,t) — A(y,t) +iB(y,t) into its real 
and imaginary parts, and get 



9A. . dB . . / dA, . dB , . \ i„ .„ 

— (y, t) + i—(y, t) = -H (-Qy-(y> t) - 'Qy-iV' *). ) + 5 ( A ' y ^ + 



(17) 
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IV. GAUSSIAN CHARACTERISTIC FUNCTION: CONSISTENT EVOLUTION 



An initial Gaussian state will keep its Gaussian form if it evolves according to a quadratic Hamiltonian dynamics. 
That means that one must expand the Hamiltonian of (jTTJ) up to order 2 in a; and y to obtain a consistent Gaussian 
evolution. Since these variables correspond to the coordinates of the trajectories supporting the state, the expansion 
will be faithful as long as the state is localized in a sufficiently small region of phase space. A coherent or even a 
squeezed state will obviously fulfill this condition in the semiclassical limit. This is the basic idea behind the following 
treatment, which can be seen as a double phase space generalization of Heller's theory of Gaussian wave packet 
evolution [10|. 

We assume that the chord function has the form 



X(y,t) =K. exp [ -a t 



■(y-Y t )-X t -h t 



1 

2fr 



(y - Y t ) ■ (r 



m t )(y-Y t ] 



(18) 



where M t and Nt are symmetric matrices. We have therefore, in the notation of the previous section, 

A(y, t) = a t -(y- Y t ) -X t + \{y-Y t )- N t (y - Y t ) 



B(y,t) = b t + -(y-Y t ) 



t(y-Y t ), 



(19) 



where Yt represents the minimum of B, or the maximum of the modulus of the wave packet. 

It is instructive to compare the above expressions with the position representation of the familiar wave packets, 
corresponding to (linearly) squeezed and rotated coherent states in [l(| [Hj]. Considering the analogy between the 
underlying double phase space coordinates (y, x) with the familiar phase space variables (g,p), we identify the matrix 
M t as describing the overall squeezing, i.e. its eigenvalues describe the compression in the chord space, y, (compensated 
by stretching in x), or vice versa. On the other hand, the matrix Nj accounts for the rotation in double phase space 
that tilts this 4-dimensional gaussian. 

Since Wj(Yu t) = 0, it is natural to expand the complex equations (JTTJ) around Y t , which leads to a separation of 
its real and imaginary parts. By identifying the terms (y — Y t ) n with n — 0, 1, 2, one obtains 



Nt 



at 



m 

dx 

dn 



dy 

■A 



N t (M t )- 1 Br t 



I —I 
dx 2 1 ' * dx 2 

d 2 n^ d 2 n J 

t — r N t -N t - 



1. 



dx 2 
-Y t -X t - 

-Y t ■ TOY* 



dx 2 



H 



d 2 n 

dydx 

d 2 n 

7T7T 1 
dydx 



d 2 H _ d 2 H 

dxdy dy 2 
d 2 H 

* dxdy 



(20) 



where Ti. and all its derivative are implicitly taken at the point (X t , Y t , t) , and we used the notations 

d 2 H d 2 H 



d 2 n 

dxdy 



dpdy p dpdy q 
d 2 H d 2 H 
dqdy p dqdy q 



= XX and 



(21) 



and 



d 2 n 

dydx 



d 2 H 
dpdy p 
d 2 H 



d 2 H 
dqdy p 
d 2 H 



dpdy q dqdy q 



= ( ^ \ 

\dxdy J 



(22) 



where T means the transpose of a matrix or a vector. 

We end up with a consistent system of ordinary differential equations. The first four equations are coupled, but 
the last two ones are actually trivial once the other are solved. Notice that these equations preserve the symmetry of 
M t and N t so we have ommited the transposition symbols that otherwise would appear in the expansion. 

In the absence of environment, that is when D = 0, the first two equations coincide with the double phase space 
trajectories of a classical dynamics, corresponding to a Liouville propagation of the chord function. The environment, 
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expressed as an exponential damping centred on y = 0, then induces a shift of the maximum of the unitary evolved 
Gaussian which, otherwise, follows those trajectories. One should note that in the latter case, the semiclassical theory 
for Gaussian evolution becomes identical to the familiar unitary theory, albeit in an enlarged phase space. Thus, the 
pair of matrices that determine the squeezing and its direction interact because of the underlying classical motion. One 
can verify that equations (|20[) are then consistent with the Linearized Green's function and wavepacket propagation 
of (Iol | , as it is partially shown in appendix [B] 

The qualitatively new feature of the present theory is that, unless D = 0, the overall amplitude of the gaussian is 
damped by decoherence, if the Gaussian is not centred on the origin. It is visible from the formal expression of M. t 
solution of (l20l) : 

rt 

JI T BQ_ T dr Q t , (23) 



(24) 



where Q t would be a time dependent 2x2 matrix solution of 

d 2 H d 2 H, 



Q t = oJ 

This leads to 



\dydx dx 2 



y-M t y= I (A ■ Q t _ T y) 2 dr, (25) 



n 



which is strictly growing. A Gaussian centred at the origin is merely squeezed by the Lindbladian term, leading to 
diffusive broadening of its Fourier transform, the Wigner function, as described in [l3l. [l5|. 

It is also interesting to go back to see how these equations read in terms of the classical hamiltonian H: 



OH ( \ \ 

K = JJ^ (X+,t) - 7 Xt + f N, + - J) M- J BJ (X+ - AT) 

X; = j|| {X-,t) -iXt+Ut- ~j) M^DJ (X+ - A7) , (26) 

where Xf = X t T \W t - Once again, when there is no environment, i.e. D = and 7 = 0, the chord tips X~^ and 
X^ are just independently following the time reverse classical motion 

a tt 

A general state can always be decomposed in a pseudo-basis of gaussian states: 

Xo(€) = / D(ij) exp [ + -^j drj. (28) 

Because of the linearity of the evolution equation of \t 1 one can evolve each of these Gaussian states independently 
by using the consistent evolution presented in this section, and the superposition of these evolved Gaussians will 
be an approximate solution of (fT2|) with the correct initial condition ([28]). Although the expression of the latter is 
written as a continuous integral, one will obviously have to use a finite, approximate decomposition. About the issue 
of optimizing the number of Gaussian states of a decomposition, one may read with interest [20| and [2l| . 

V. QUADRATIC HAMILTONIAN 

Here we discuss the consistent Gaussian evolution in the case of a quadratic Hamiltonian, namely 

H(x)=x-Wx, (29) 
with some symmetric matrix EL The double phase space Hamiltonian is then 

H(x, y) = -2x ■ HJJy - 7 x ■ y. (30) 
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It is obvious that the consistent evolution will give an exact result here, since the Hamiltonian coincides with its 
quadratic expansion. However it is instructive to explicit the calculus to have an insight on the role of the different 
terms. 

As a reminder, the general solution derived in [13l ] was written in the following form 



-tt 



£= -R-Jy = -l[Rjy) ■ (34) 



H 



e~ 27 *R f MM t T + / e^'-^Rt-t'^-J-t' dt' 



(31) 



where Xo(£) is anv general (not necessarily localized) initial function, and 

R t = exp(2JJEft). (32) 

With the notation adopted in this paper, it reads 

A(y,t) = A (y )=A (e-"< t Rjy) 

B{y, t) = B (e-^Rjy) + ~ J e 2 ^'^ (a • K t T _ t ,y) 2 dt' . (33) 
Here we used the correspondence 



If the initial state is a Gaussian wave packet, it imposes the form of Ao(y, t) and Bo(y, t), and one has 
A(y,t) = a-(e-^ t Rjy-Y)-X + ^(e-^Rjy-Y)-N(e-' t Rjy-Y) 

B{y, t) = b+ \{e-^Rjy - Y) ■ M(e~^R[y - Y) + I J e 2 ^ 1 '^ (a • Rj_ t ,y) * dt' . (35) 

This solution can be matched with (|19p by setting the different parameters in the following way, 

N t = e _27 *R t NR7 



X t = e -7 R t (X + NY) — N t Y t 
Y t = (Mty 1 e-^RMY 

a t = a + X ■ Y + Iy ■ NY - X t ■ Y t - W t ■ N t Y t 

bt = b + iy • M • Y - W t ■ M t Y u (36) 

where the matrix D is set for ^AA T ^j; Yt is defined as the maximum of the Gaussian, ^(Y t ,t) = 0, which is not 

the Hamiltonian evolution of the initial maximum Y of the Gaussian; and X t is defined as —^(Y tl t), that is, the 
argument of the order of the expansion of H in (TIT)) . 

On the other hand we can explicit the consistent Gaussian evolution given by equations (|20[) : 



x t 


= (2JH- 


7 ) Xt + NtiMt^BYt 


Y t 


= (2HJJ + 


7 ) Y t - (M t )- x By t 


N* 


= 2JJHN t 


- 2N t HJ - 2 7 N t 


Mt 


= 2JJHM t 


- 2M t HJJ - 2 7 M t + D 


dt 


= X t -(W 

= Iy t -i 

2 ' 


hT^Yt 


bt 





(37) 

One can then check that (f36|) is solution of (f37|) . 

The equations for the evolution of the coefficients of the Gaussian is far less intuitive than the simple interpretation 
of the solution (f3T|) . The latter is indeed a product of two Gaussians, each of which is fairly simple. The product of 
two Gaussians is itself a Gaussian, but its maximum does not obey any simple differential equation. Instead of fllSp . 
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one could consider using an a priori solution structured as a product of two Gaussians, and perform the expansion of 
the Hamiltonian Ti of (fl7| , not around Y = Y t , but around both Y = and Y = e~ 7t RjY. However, none of these 

points correspond to a zero of §pg(Y, t) to justify the expansion around Ti (^— §pg(Y, t), tj , which permits to separate 

the real from the imaginary part of (|17p in the Gaussian approach. This is not a problem in the quadratic case where 
the expansion is exact whatever the origin: 

d A ft R ft A ftH ? 

— (y, t) + i^(y, t) =-^-(y,t). (2HJJ + 7 ) y - i — (y, t) ■ (2HJ + 7 ) y + -(A ■ y) 2 . (38) 

In the non-quadratic case, this exact separation does not occur any more, and one must then relie on an awkward 
dynamics for Y t and X t . However a classical dynamics can still be obtained by performing a complex WKB analysis, 
which will be addressed in 



VI. CONCLUSION 



We have derived an approximate Gaussian solution of the Lindblad equation, based on a second order expansion 
of the Hamilton Jacobi equation of the phase of the chord function. The trajectory of the maximum of the Gaussian 
does not follow an intuitive trajectory, as it is a compromise between the unitary evolution of the Gaussian, and the 
damping induced by the Lindbladian term. The latter quickly shrinks the relevant part of the chord function to a 
region close to y = on the y/h scale, which can be interpreted as the region containing all the classical information 
about the state. Thus the maximum of the Gaussian chord function ends up being at y = 0. This result is in 
agreement with established results in the case of a quadratic Hamiltonian and linear Lindblad operators, see for 
instance [l3j], as well as fulfilling our handwaving intuition that "decoherence drives the system back into classical 
dynamics" . However, with a nonquadratic Hamiltonian, the rate of this environmentally induced exponential damping 
acquires a phase space dependence. 

Our theory generalizes the standard semiclassical approximations for the unitary evolution of wave packets [13 ] in 
a nontrivial way: in the chord representation adopted here, a single wave packet is always centred on the origin (its 
position is represented by an overall phase). Thus, by allowing gaussians away from the origin, like in the typical 
case of a "Shrodinger cat" state recalled in appendix [C] the approximate gaussian evolution described by the theory 
includes the interferences between pairs of wave packets, their motion and their loss of coherence. 



APPENDIX A: ASYMPTOTIC EXPANSION OF THE HAMILTONIAN OPERATOR 



One can verify easily that 



II) (31 • vx],xs 



dSyfdSy 1 [ m(m - 1) d 2 S fdS 
dp J \dq J A 

d 2 s fds^™- 1 



2q /aq\ m-2 /gs 



dpdq \ dp 



dq 



2 dp 2 \ dp 

n(n-l) d 2 S fdS 



2 dq 2 \ dp 



8S 
dq 



+c 4 



exp XS. 



From there one can generalize and write 



(Al) 



n 1 --iW V ) 6XP fa ' = 



H I —z-,v 
dy 



2% X dx 2 



dS \ d 2 S, 



exp , -S 



(A2) 



APPENDIX B: CONSISTENCE WITH HELLER'S EQUATIONS 

In the paragraph 2.3. Linearized Green's function and wavepacket propagation of [13], Heller starts with a Gaussian 
wave function with 



Ml, t) = exp {i/h[A(q - Q) 2 + P(q -Q) + s t }}, 



(Bl) 
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where A — a + ib is a time dependent complex number, P a real time dependent number and we have sticked to the 
1 dimensionnal case. Then he proceeds with 



El 

2ni 



Q = 
P = 
A = 

Pz = 
Z = 

V(Q), and V" for 



dH 

dp 
dH 



El 
2z 

— V Z 
1 

— Pz, 

m 



0W)- 



(B2) 



This defines the consistent Gaussian dynamics of 



where H stands for H(P, Q) 
(|B1I) according to Heller. 

In this appendix we show that the evolution of the chord representation Xhif) of the state \i/jh,(t))(iph(t)\, as defined 
above, also follows equations (|2"0)) in the unitary case D = 0. This shows in a simple case that our consistent Gaussian 
evolution is consistent with the one of Heller in the unitary case. Both theories also agree for a general state, although 
the derivation, carrying much heavier expressions, is not written out in this appendix. 

From (|B1[) and ([5]) one has 



Xh(t) = exp- 



-PVp 



with 



QVq ' 








a 2 


a 


+ T 


¥ 


a 




2l> 


-Lb 



-V ■ My 



(B3) 



(B4) 



Notice already that the first two lines of (|B2[) give the second line of ([20)1 with X t = (P,Q). On the other hand 
Y t = and N< = in this case. Moreover, the equation for A in (|B2[) gives 



/ ab- — _ + - a 2 -6 2 

m b \ 2 m 



2b 2 



b) m ab -2b {— 
— ab 



-{a 2 -b 2 - 



JL± ab \ 



2b 2 



1 4 

Ab 2 m 



ab 



(B5) 



that is 



—V' 



-2V"I 



"■qq 



<Lqq 

1 R 



-V"l 



It is equivalent to the fourth line of (|20|1 . with 



d 2 n 

dydx 



r i><i 



-V"{Q) 




(B6) 



(B7) 



We do not take into account St which corresponds to the prefactor of \, which would be given by the next order in h 
in our expansion. 

To check the general case, with Y t ^ and Nt ^ 0, one must write the chord representation x'ab of the cross product 
\ipa(t))(ipfr(t)\ of two different Gaussian states which follow equations (|B2[) . The expression of the time derivative of 
that chord function can be shown to be consistent with (|2"01) in the same way as what is done above, but quite heavier. 



APPENDIX C: "SCHRODINGER CAT" STATES 



The generalized "cat state" , that is a coherent superposition of two squeezed states, is an important example for 
two reasons. First it is Gaussian, and therefore is adapted to the treatment of this article; and secondly it gives a 
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quite transparent physical interpretation of the chord function. We define the two ingredients of such a state: two 
squeezed states 

Mq) = ,9M/4 eX P 



(^2)1/4 *y 2uJ 2 h h 
i . / {q-Qbf , iPi 



Mq) = p^ exp (r^I?r + -?r> (C1) 

centred respectively on X a = (P a ,Q a ) and X b = (Pb,Qb)- The density operator of the corresponding "cat state", 
\4>a) + IV'b); has two diagonal terms and two off-diagonal terms: 

p= \\i>a)m + Im^bi + \mm + m 

We are interested here in the sum of the pair of non-local terms, \ip a }(ipb\ and complex conjugate. Notice that the two 
other diagonal terms, which have a classical interpretation, can be retrieved anyway by setting a = b. In the Weyl 
representation, this non-local term can be written 



1 2u a u b ( (q-Qf+^l(p-Pf+l^l-Lol)(p-P){q-Q) 



(-V ■ n 

h 



+ l (-Y-x-PY p )\, (C3) 



where X = (P, Q) = — *± — - and JY = X b — X a . In the chord representation, one has 



1 / l^b ( (y P - Y p f + wfa*(y q - Y q f + »K 2 - o;g)(y p - y p )(y - Y q ) 

Xab(V) - TyZTfc \ / , .2 , , .2 CX P ~ """" """ """"""" 



i 



Then one has 



and 



+-(-X.y + Y 9 Q)). (C4) 



Mv) = \ i ^tr\ { yv - n»)(y, - n) - * ■ v + W (C5) 

ZuJ a+ UJ b 



One may have noticed that A (y) is not an odd function, nor B (y) an even function, but one has to bear in mind 
that we looked only at Xab{y), and that the complete chord function contains also Xba{y), which re-establishes the 
overall symmetry, x{y) — x{~y)* f° r Hermitian operators. 
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